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^ ABSTRACT 



The problem of recovering geometric properties of a domain from the trace of the heat 
kernel for an initial-boundary value problem arises in NMR microscopy and other appli- 
O ! cations. It is similar to the problem of "hearing the shape of a drum" , for which a Poisson 
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type summation formula relates geometric properties of the domain to the eigenvalues of 



O ' the Dirichlet or Neumann problem for the Laplace equation. It is well known that the 
^ ■ area, circumference, and the number of holes in a planar domain can be recovered from 
. the short time asymptotics of the solution of the initial-boundary value problem for the 
P_i! heat equation. It is also known that the length spectrum of closed billiard ball trajec- 
^ I tories in the domain can be recovered from the eigenvalues or from the solution of the 
+^ ' wave equation. This spectrum can also be recovered from the heat kernel for a compact 
g ■ manifold without boundary. We show that for a planar domain with boundary, the length 
spectrum can be recovered directly from the short time expansion of the trace of the heat 
kernel. The results can be extended to higher dimensions in a straightforward manner. 



^Partially supported by a research grant from the Foundation for Basic Research administered by the 
Israel Academy of Science and by a research grant from the US- Israel Binational Science Foundation. 



1. Introduction 



The problem of recovering geometric properties of a domain from NMR measure- 
ments arises in oil explorations and in non- invasive microscopy of cell structure pp. In 
these measurements the trace of the heat kernel for the initial value problem with reflect- 
ing (Neumann) boundary conditions is measured directly. The problem is analogous to 
"hearing the shape of a drum" , where the solution of the wave equation in the domain is 
measured directly (it is "heard"). 

The problem of recovering geometrical properties of a domain from the eigenvalues of 
the Dirichlet or Neumann problem for the Laplace equation in a domain has attracted 
much attention in the literature (see [21- [Z] for some history and early results; for more 
recent work see [H], an references therein). 

The mathematical statement of the problem is as follows. Green's function for the heat 
equation in a smooth planar domain Q, with homogeneous Dirichlet boundary conditions, 
satisfies 

^^^^'^"^'^^ = DAyG{y,x,t) foTy,xen,t>0 (1.1) 
G{y,x,0) = S{y-x) (1.2) 



G{y,x,t) = for y E dn, X en, t > 0. (1.3) 

The function G{x, x, t) dx is the probability of return to x dx at time t of a free Brownian 
particle that starts at the point x at time t = and diffuses in Vt with diffusion coefficient 
1, with absorption at the boundary dVL. If it is reflected at dVt, rather than absorbed, the 
Dirichlet boundary condition ()1.3p is replaced with the Neumann condition jTUI 

"^^i^;"^'^^ =0 ioiyedn,xen,t>Q, (i.4) 

du[y) 

where ^{y) is the unit outer normal at the boundary point y. The trace of the heat kernel 
is deflned as 

P{t) = f G(x,x,t)dx (1.5) 
Jn 

and can be represented by the Dirichlet series 

oo 

Pit) = Y.e-'"\ (1.6) 



n=l 



where are the eigenvalues of Laplace equation with the Dirichlet or Neumann boundary 
conditions p.3|) or p.4|l . respectively. 
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It has been shown by Kac [2] that for a domain Q with smooth boundary dQ, the 
leading terms in the expansion of P(t) in powers of \/t are 

W)~§-|^ + i(l-'-) + 0(v^). for^^O. 

where \ denotes the area of Q, \dQ\ denotes the arc-length of dQ, and r is the number 
of holes in Q. The full short time asymptotic power series expansion of P{t) in the form 

oo 
n=0 

can be deduced from the large s expansion of the Laplace transform 

g{s) = exp{-.2t} [P{t) - ^) dt, [a, = M j . 

in inverse powers of s. Such an expansion was given by Stewartson and Waechter [Hj in 
the form 



°° c 



where 



n=l 



ar,r 



The constants c„ are computable functionals of the curvature of the boundary. The full 
expansion is denoted 

If the boundary is not smooth, but has cusps and corners, the expansion contains a term 
of the order t~'^, where u is a number between and 1/2. 

The Stewartson- Waechter expansion was used in |8J to deduce further geometric prop- 
erties of fl by extending g{s) into the complex plane. Examples were given in ^ of the 
resurgence of the length spectrum of closed billiard ball trajectories in the domain. 

The full length spectrum of closed geodesies on a compact Riemannian manifold with- 
out boundary Q appeared in the short time asymptotic expansion given in jB], 

1 OD 

P(t)r~. \2PniVi)e-'-/\ fort^O, (1.8) 

where 6n are the lengths of closed geodesies on Q and Pn{x) are power series in x. 

In this paper, we construct an expansion of the form ()1.8|1 for the trace of the heat 
kernel for the initial-boundary value problem (jl.lll - fjl.Hj) or ()1.4j) in a smooth bounded 
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domain in M^. The results can be generalized to higher dimensions in a straightforward 
manner. 

The point of departure for our analysis is the observation that transcendentally small 
terms are not included in the expansion ()1.7j) . These terms have been neglected in [H] and 
|Sj even in the case of a circular domain, where the Laplace transform of G{y^ x, t) can be 
expressed explicitly in terms of modified Bessel functions. In |S| this Laplace transform 
is expanded in inverse powers of s and the coefficients are evaluated asymptotically for 
large n. 

A generalization of the asymptotic expansion "beyond all orders" p.8|l has the form 
P{t)r^Pswit) + ^,J2Pn{Vi)e-''-/\ fort^O, (1.9) 



n=l 



where 6n, ordered by magnitude, are constants to be determined, and Pn{x) are power 
series in x. Transcendentally small terms may be, in fact, quite large and make a finite 
contribution to the expansion ()1.9p jTT]. 

To recover the geometrical information from the expansion ()1.9jl . given the (measured) 
function P{t), we note that 



\n\ = lim47rtP(t), \dn\ = -liniSVTrt 



P{t)- 



(1-10) 



and so on. This way the entire expansion (jl.7p can be determined. 

Once the coefficients of the expansion ()1.7p have been determined, the exponent of the 
dominant term of the transcendentally small part, 6i, is found as 

5i = -limtlog [P(t)-PswW 

Proceeding this way, we can recover the entire expansion (|1.9p if P{t) is known (e.g., from 
measurements). 

In this paper, we use the "ray method", as developed in to construct a short 

time asymptotic expansion of the heat kernel. We use it to expand the trace asymptot- 
ically beyond all orders (the so called "hyperasymptotic expansion") and show that the 
exponents 6i are the squares of half the lengths of the periodic orbits in the domain. The 
exponentially small terms in the expansion ()1.9p are due to rays reflected in the bound- 
ary, much like in the geometric theory of diffraction jn]-|IH]- This recovers the length 
spectrum of closed billiard ball trajectories in the domain. In particular, the smallest 
exponent 6i is the width of the narrowest bottleneck in the domain. 

2. The one-dimensional case 

The solution of the heat equation in an interval can be constructed by the method of 
images. Speciflcally, the Green function of the problem satisfles 

dG{y,x,t) d'^G{y,x,t) , „^ ^ 

di " dy^ for < X,?/ < a, t > (2.1) 
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G{y,x,0) 



— 5{y — x) ior < x,y < a 



(2.2) 



G'(0,a;,t) 



G{a, x,t) = for < X < a, i > 0, = 0, 1. (2.3) 



The method of images gives the representation 

{y — X + 2naY \ 



G{y,x,t) 
1 



2V7rt 



E 



exp 



At 



I)'' exp 



{y + X + 2nay 
At 



(2.4) 

, (A; = 0,1). 



Note that if the infinite series is truncated after a finite number of terms, the boundary 
conditions are satisfied only in an asymptotic sense as i — > 0. That is, the boundary values 
of the truncated solution decay exponentially fast in t~^ as i — > and the exponential 
rate increases together with the number of retained terms. 
The trace is given by 



G{x, X, t) dx 



1 



2VTitJo 

1 oo 

—— T 



E 



\na\ 



exp 



n=— oo 
oo 



\na] 



a exp 



| + (_l)^exp| 

+ (-1)^= Texp 
Jo 



[x + na) I 

{x + naY 



dx 

dx 



2Vnt 



E 



Ina) 



+ 



-11 



a ^ (-1 



2V7rt 2 2V7rt ^0 



E^xp 



ina] 



On the other hand. 



r-a oo 

/ G{x,x,t)dx =J2(^~^"^^ 

•'^ n=l 



, (A; = 0,1). (2.5) 



(2.6) 



where {An} are the eigenvalues of the homogeneous Dirichlet or Neumann problem for 
the operator d'^/dx'^ in the interval [0, a]. Thus 



E' 

n=l 



2Vnt 



+ 



+ 



a 



2V7rt 



^exp 



(na)^ 



, {k^O, 1] 



(2.7) 



If instead of a single interval of length a, we consider the heat equation in a set Q 
consisting of K disjoint intervals of lengths Ij, {j = 1, .., K), respectively, the resulting 
expansion is 



Ee 

n=l 



= + (-1)'^ + E E exp 



, (A; = 0,1). (2.8) 



The numerator in the first term on the right hand side of eq. ()2.8p can be interpreted 
as the "area" of fi, so we denote it Y,f=i^j = 1^1- The number 2K is the number of 
boundary points of Q, which can be interpreted as the "circumference" of the boundary, 
so we denote it \dQ\ = 2K. The exponents in the sum on the right hand side of eg. ()2.8j) 
can be interpreted as the "widths" of the components of Q. Clearly, for small t, the term 
containing the smallest width, r = mini<j<i^ /j, will dominate the sum. Thus we can 
rewrite eq. ()2.8|) as 

E e = 777= - V + ^ -T + E ;H=7 E —T^ ' (2-9) 



n=i 2v%t 4 Vvrt [ t } 1% 2Vvrf ^0 [ t J 

where m is the number of the shortest intervals in Q. 

Equation ()2.9|1 can be viewed as the short time asymptotic expansion of the sum on 
the left hand side of the equation. The algebraic part of the expansion consists of the first 
two terms and all other terms are transcendentally small. The geometric information in 
the various terms of the expansion consists of the "area" of Q and the "circumference" 
\dfl\ in the algebraic part of the expansion. The transcendental part of the expansion is 
dominated by the term containing the smallest "width" of the domain, r. 

The geometric information about Q contained in the algebraic part is the information 
given in the "Can one hear the shape of a drum" expansions |2], j^]- The geometric 
information contained in the transcendentally small terms in ()2.9p can be understood as 
follows. The terms nlj in the exponents are the lengths of closed trajectories of billiard 
balls in Q, or the lengths of closed rays reflected at the boundaries, as in "Sj. 

The representation ()2.4|) can be constructed as a short time approximation to the 
solution of the heat equation ()2.1|) - ()2.3p by the ray method [12]. In this method the 
solution is constructed in the form 

G{y, X, t) = ^ Z^y, x)e~'l\ (2.10) 

n=0 

Substituting the expansion ()2.10|) into the heat equation ()2.H) and ordering terms by 
orders of magnitude for small t, we obtain at the leading order the ray equation, also 
called the eikonal equation, 



dS{y,x" 



dy 

and at the next orders the transport equations 



(2.11) 



dS{y,x)dZ^{y,x) (d^S{y,x) 2n 

+ '^n[y, X) h 



dy dy \ dy^ S{y,x)^ 

2 d^Zr,.i{y,x) 



S{y,x) dy"^ 



n = 0,l,.... (2.12) 
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Denoting 

dS{y,x) 

we write the equations of the characteristics, or rays of the eikonal equation ()2.11|) as JH] 

with the initial conditions 

?/(0, x) = X, p(0) = ±1, 5(0) = 0. 

The condition S'(O) = is imphed by the initial condition G{x,y,Qi) = S{x — y). The 
solutions are given by 

y(r, x) = X + 2pr, p(r) = ±1, 5(r) = 2r = ±(y - x). (2.14) 

Thus S{y,x) is the length of the ray from y to x. We denote this solution by So{y,x). 
It is easy to see that the solution of the transport equations corresponding to So{y,x) is 
given by Zo(|/,x) = const, and Zn{y,x) = for all n > 1. The initial condition ()2.2j] 
implies that 

2V7r 

Combined in eq. ()2.10|) this solution gives Green's function for the heat equation on the 
entire line, 

which is the positive term corresponding to n = in the expansion ()2.4|) . 

The ray from x to ?/ is not the only one emanating from x. There are rays emanating 
from X that end at y after reflection in the boundary. Thus the ray from x that reaches y 
after it is reflected at the boundary has length y + x. Therefore there is another solution 
of the eikonal equation, Si{y,x), which is the length of the reflected ray, given by 

Siiy,x) = y + x. 

The ray from x that reaches y after it is reflected at the boundary a has length 2a — x — y. 
The ray from x to 0, then to a, and then to y has length 2a + x — y. Thus the lengths of all 
rays that reach y from x after any number of reflections in the boundary generate solutions 
of the eikonal equation, which are the lengths of the rays, which in turn generate solutions 
of the heat equation. We denote them by Sk{y, x) with some ordering. The corresponding 
solutions of the transport equation are 
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where Ck are constant. They are chosen so that the sum of all the ray solutions, 



satisfies the boundary conditions ()2.3|) . Note that for all A; 7^ 

Gk{y, X, t) — as t — s> 0. 
This construction recovers the solution ()2.4|) . 

3. The ray method for short time asymptotics of Green's function 

The ray method consists in the construction of Green's function G{y, x, t) in the 
asymptotic form 



n-l 



(3.1) 



n=0 



The function S{y, x) is the solution of the eikonal equation 



^ySiy,x] 



and the functions Zn{y, x) solve the transport equations 



(3.2) 



2VyS{y, x) ■ VyZ,,{y, x) + Z.^{y, x) 
2 



AyS{y, x) + 



2n - 1 



S{y,x) 



AyZn-i{y,x), forn = 0, 1,2, . . .. 



(3.3) 



The eikonal equation (j3.2j) is solved by the method of characteristics ^3] ■ The char- 
acteristics, called rays, satisfy the differential equations 

'-^ = 2VySiyir,x),x), dVyS{yir,x),x) ^ dS{y{r,x),x) 



dr 



dr 



The initial condition ()1.2|) implies that the rays emanate from the point x. Thus we 
choose the initial conditions 

y{0,x) = x, VySiy{0,x),x) = u, S{y{0,x),x) = 0, (3.5) 

where 1/ is a constant vector of unit length. The solution is given by 

y{T,x) = X + 2iyT, S{y,x) = \y - x\ = 2t, VyS{y,x) = u. (3.6) 
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The pair (r, u) determines uniquely the point y = y{T, x) and the value of S{y, x) at the 
point. The parameter r is half the distance from to cc or half the length of the ray from 
X to y. The vector u is the unit vector in the direction from x to y. 

The function Zo{y,x) is easily seen to be a constant, l/47r, and Zn{y,x) = for all 
n > 0. This construction recovers the solution of the heat equation in the entire plane 
and disregards the boundary dQ, because in the plane every point can be seen from every 
other point by a straight ray. Note that to calculate the function P{t) in eq. ()1.5p only the 
values of S{x, x) and Zq{x, x) are needed. Thus S{x, x) = and the first approximation 
to G{x, X, t) is 



hence the first approximation to P{t) is 

Poit) 



47d' 



There is another solution of the eikonal equation ()3.2|) constructed along rays that 
emanate from x, but reach y after they are reflected in dfl The law of reflection is 
determined from the boundary conditions. Dirichlet and Neumann boundary conditions 
imply that the angle of incidence equals that of reflection ^2] • Similarly, there are solu- 
tions of the eikonal equation that are the lengths of rays that emanate from x and reach 
y after any number of reflections in dfl. We denote these solutions Sk{y,x) with some 
ordering. Thus the full ray expansion of Green's function has the form 



G{y,x,t) ~ ^e-^^(?/'^)/4^Z,(t/,a;,t), (3.7) 
fe=i 

where 

oo 

Zk{y,x,t) = ^ Z„_fe(y,a;)r-\ 

n=0 

As above, each one of the series 

e-Sl{y,x)/At2^^y^^^f^ 

is called a ray solution of the diffusion equation. The boundary values of Zk{y,x,t) 
are chosen so that G{y,x,t) in eq. (j3.7l ) satisfies the imposed boundary condition. In 
particular, the values of Sk{x, x) are the lengths of all rays that emanate from x and are 
reflected from the boundary back to x. Note that sums of ray solutions satisfy boundary 
conditions only at certain points. 

To fix the ideas, we consider first simply connected domains. We denote 

Soiy,x) = \x-y\ 
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and 

G„(i,,.=.t) = ^e-^o(f.-)/«. 

We consider first solutions corresponding to rays that are reflected only once at the bound- 
ary, and in particular, rays that are reflected back from the boundary to the points of 
their origin. Such rays hit the boundary at right angles (see Fig. 1 and If there 

is only one minimal eikonal Si{x, x) > 0, we say that a; is a regular point of fl. If there 
is more than one minimal eikonal Si{x,x), we say that x is a critical point of fl. We 
denote by F the locus of critical points in Q. The eikonal Si{y,x) is the length of the 
shortest ray from x to y with one reflection in the boundary such that the ray from x to 
the boundary does not intersect F. For x = y the eikonal 6*1(0;, a;) is twice the distance 
of X to the boundary. We denote by x' the orthogonal projection of x on the boundary 
along the shortest normal from x to the boundary. When y = x' 

Si{x',x) = So{x',x) = \x-x'\. (3.8) 

The function 

G,iy,x,t)=e-S"(y^^y^tz^ {y,x,t) 
has to be chosen so that Go{x', x, t) —Gi{x', x, t) = 0. In view of ()3.8p . we have to choose 

Z^ (aj>,t) = -i-. 

When y" is the other boundary point on the normal from x' to x, we have 

Goiy", X, t) - G,{y\ x, t) = ^^e-l^-?/"!^* - ,-i\x'-xiMy"-x'if/tz^ (^.^ ^ (3 9) 

Next, we consider in f2 — F the minimal among the remaining eikonals Sk{x,x) > 
Si{x, x) and denote it S2{x, x). This eikonal is twice the length of a ray that emanates 
from X, intersects F once, and intersects the boundary dfl at right angles at a point, 
denoted x". The eikonal 6*2 (y, x) is the length of the ray from x to y with one reflection 
in the boundary such that the ray from x to the boundary intersects F once. When 

y = x" 

Siix", x) = So{x", x) = \x- x"\ . (3.10) 

When y' is the other boundary point on the normal that emanates from x" (see Fig. 2), 
we have 

O f I ^\ I "III' „"l 

b2[y ■, x) = \x - X \ + \y - X I . 

In general x' 7^ y' and x" 7^ y" . However, if the ray is a 2-periodic orbit (that hits the 
boundary at only 2 points), x' = y' and x" = y" so that 

S2iy",x) = Soiy",x) = \x-y"\ 

and 

S2{y',x) = \x-x"\ + \y" -x'\. 
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Figure 1. The locus of critical points, T , is the segment AB. The first eikonal is 
Si{y, x) = \x — c\ + \c — y\. It is defined as the shortest reflected ray from x to y, such 
that X — c does not intersect T . For x = y the diagonal values are Si{x,x) = 2\x — x'\. 
The diagonal values of the second eikonal are "\. The vectors x — x' 

and X — x" are orthogonal to the boundary. For x\ the two eikonals are equal. 

Since 

I „ "1^1 'III" '1^1 III" 'I 

\x — y \ < \x — X \ + \y — X \ < \x — X \ + \y — x \ 

for all regular points x, the order of magnitude of the boundary error ()3.9|) decreases if 
we use the approximation 

Go{y, X, t) ~ Go{y, x, t) - Gi{y, x, t) - Gad/, t) (3.11) 

with 

Z2{y",x,t) = Z,{y",x,t) = Zo{t). 
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Figure 2. The second eikonal 6*2(2/, x) = \x — d\ + \d — y\. It is defined as the shortest 
reflected ray such that x — d intersects T . The eikonals S^lx, x) and 5*4(0;, x) are ordered 
according to magnitude. 

4. The trace 

To find the short time asymptotics of the Dirichlet series as given in eq. (|1.5|) . 



we use the ray expansion ()3.7|) for the evaluation of the integral. We retain in the result- 
ing expansion only terms that are transcendentally small, since all algebraic terms are 
contained in the expansion ()1.7|1 . 

4.1 Simply connected domains 

We note that according to Sard's theorem, F is a set of measure zero and that all 
points in the domain f2 — F are regular. For any point x E Q, we denote by ri{x) its 
distance to the boundary and note that Si{x,x) = 2ri{x). We also denote by Si{x) 
the arclength at the boundary point x' (the orthogonal projection of x on dQ along the 
shortest normal from x to dQ), measured from a boundary point where the arclength is 
set to (see Figure 3). 
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Figure 3. The arclength Si{x) is measured from the point E. Both transformations 
X {T^ii^), ^1(3^)) o.iT'd X — > {r2{x), Si{x)) are one to one mappings o/Q — F. The images 
are given in Figure 4- 

It follows that the change of variables in O — F , given by 

{ri{x),si{x)), (4.1) 

is a one-to-one mapping of ^2 — F onto a strip < ri < ri (si) , < Si < L, where ri(si) 
is the distance from the boundary point corresponding to arclength si to F. 
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Figure 4- The domain Q is the ellipse ^ + ttt < 1- The domain enclosed between the 

Si-axis and the lower curve is the image of the ellipse under the transformation \4.1^ and 
the domain enclosed between the upper and the lower curves is its image under 

We evaluate the integral over Q separately for each summand k in the expansion ()3.7p . 
In this notation, we can write 



/ Gi{x,x,t) dx 
Jn 



(4.2) 



Si{x,x] 



/4t ^ rri{si) 2 /, 

22 Zn,iix,x)t''^^ dx = ds e-''i/Vi(ri,si)Zi(ri,si,t) rfri 

n=o -^0 ^0 



where Ji(ri, si) is the Jacobian of the transformation and 



Zi(ri,si,t) = Y.Zn,iix,x)t''-\ 

n=0 

Note that the Jacobian vanishes neither inside f2 — F nor at ri = 0, because the transfor- 
mation is one-to-one in f2 — F, however, it does on F. 

We set S2{x,x) = 2r2{x) and use it as a coordinate. We use si{x) as the other 
coordinate of the point a; G — F. Note that while r2{x) is the length of the longer 
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normal from x to dQ (the one that intersects F), the other coordinate is the arclength 
corresponding to the shorter normal from x to dQ (the one that does not intersect F). 
The transformation 

x^{r2ix),siix)) (4.3) 

maps Q — T onto the strip r (si) < r2 < /(si), < Si < L, where /(si) is the length of the 
segment of the normal that starts at the boundary point ri = 0, Si and ends at its other 
intersection point with the boundary. This mapping is one-to-one as well. It follows that 



/ G2(x,x,t)dx 
Jn 

S2ix,x) 



(4.4) 



/4t 



In 
where 



J2Zn,2{x,x)r-Ux 
n=0 



L Msi) ^ 2,, 

dsi / e '"2/ J2ir2, Si)Z2{r2, si,t) dr2, 

Jr{si) 



Z2{r2,s^,t) = Y.Z^,2{x,x)r-\ 

n=0 

Note that for a; on F both transformations ()4.1|) and ()4.3|) are identical and 

J2ir2,Si)Z2ir2,Si,t) = Ji(ri,Si)Zi(ri,Si,t). 
It follows that the two equations ()4.2|) and ()4.4|) combine together to give 



/ \Gi{x,x,t) + G2{x,x,t)] dx 
Jn 



(4.5) 



r /'^'^ e-"'/* J(r, s)Z{r, s, t) dr ds, 
Jo Jo 



where s = Si, r = ri, J = Ji, and Z = Zi for < r < ri(si), and s = Si, r = r2, J = J2, 
and Z = Z2 for r2(si) < r < l{si). Thus the domain of integration of the function 
g-r /tj(^^^ s)Z(r, s, t) in eg. ()4.5|1 is the domain enclosed by the Si-axis and the upper curve 
in Figure 4. Now, for t -C 1, we write the inner integral on the right hand side of eg. ()4.5j) 

as 



e-' /Kj{r,s)Z{r,s,t) 



dr = ^erf (^-^^ J(0, s)Z{0, s, t) (l + O [Vt 



exp 



/(.) 



j(o,s)z(o,s,t) {i + o(yt 
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Recall that J(0, s)Z{0, s, t) ^ 0. Only the exponentially small terms have to be considered, 
because the algebraic terms are included in the SW expansion. Thus 

e-"'/* J(r, s)Z{r, s, t) dr ds - \/y (0, s)Z(0, s, t) (l + O (Vt)) ds = 
P{s) \ J{0,s)Z{0,s,t) 



exp 



-O (t) ds for t < 1. 



Evaluating the last integral by the Laplace method, we find that each point Sj that is an 
extremum point of l{s) contributes and exponential term of the form 



exp 



l\s,)\ J{0,s,)Z{0,s„t) 



o in . 



(4.6) 



The expression (|4.6jl means that some of the 5„-s in the expansion eq. (|1.9|) are the extremal 
values l{si) and their multiples. These are half the lengths of the 2-periodic orbits of a 
billiard ball in Q (see Figure 5). The 2-periodic orbits of the ellipse are the major axes, 
which correspond to the lowest and highest points of the top curve in Figure 4. There are 
other exponents as well, as discussed below. 




Figure 5. The rays emanating from the boundary points Si and S2 are orthogonal the 
boundary at both ends. They are 2-periodic orbits. 



The pre-exponential terms in the expression ()4.6|) influence the factors Pn{\/i) in 
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eq. ()1.9|) . For example, if l'{si) = 0, l"{si) ^ 0, then v = 3/2. If the boundary is flatter, 
then 1 < < 3/2. 

In addition to the 2-periodic orbits, there are ray solutions corresponding to rays from 
X to y that are reflected any number of times in the boundary. There are eikonals from x 
to 2/ in f2 with N — 1 different vertices on the boundary, which have vertices on dQ if 
X = y and x G dfl (this is a periodic orbit with — 1 reflections). Among these periodic 
orbits there are eikonals Siy{x, x) with extremal length, denoted S^j, (j = 1, • • •)■ 
points X & Q on a 2-periodic orbit the eikonal Sj\f{x,x), which now has A^ — 1 vertices 
on the boundary, may reduce to the 2-periodic orbit with A^ reflections. Therefore the 
change of variables x —>■ {Sn{x, x), s{x)) will map the domain into a strip with extremal 
widths that are the differences between the lengths S^j and the length of a 2-periodic 
orbit with A^ reflections. It follows that the evaluation of the trace by the Laplace method 
leads to exponents which are the extremal lengths of periodic orbits with any number of 
reflections. 

For example, there is an eikonal in a circle (centered at the origin) that is the ray from 
X to y with 2 reflections in the boundary (see Figure 6). 




Figure 6. The eikonal Sz{y,x) with two reflections in the circle. 
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For a; = 7/ it is the equilateral triangle (see Figure 7) with circumference 



S{x,x) = R 



2 a; 2 + 1 + \ 8 x ^ + 1 



V 



4|a;|2 + 1 + J8\x 



+ \/4 a;P + 2 + 2A/8 a;P + 1 









R 








|x| 


P 








\ /A" 

\^ ■ a \ 






fa ix/. 




\ 







/ 



Figure 7. The eikonal S^{x,x) with two reflections, where \x\ = OC . 



The eikonal S^{x,y) reduces to a 2-periodic orbit with two reflections ii x = y = Q 
(the center of the circle). If x is on the circumference, the eikonal becomes the isosceles 
triangle with one vertex at a;. To evaluate the contribution of the corresponding ray 
solution to the trace, we use this eikonal as a coordinate that varies between 4i?, the 
length of the 2-periodic orbit with two reflections, and 3v^-R, the circumference of the 
inscribed isosceles triangle. The contribution of this integral to the exponential sum in 
eq. ()1.9|) contains exponents that are both lengths. 

Similarly, the 2-periodic orbit with 3 reflections has length 6i? while the periodic orbit 
with 3 reflections at 3 different points has length A\/2R < 6R. 

4.2 Multiply connected domains 

Once again, we consider flrst rays from x to y that are reflected only once in the 
boundary. For every connected component of dQ, denoted dfli (i = 1, ...,/), a point x 
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in Q is regular with respect to dQi if there is only one minimal eikonal Si^i{x, x) > with 
one reflection at dQi. We denote by Fj the locus of the irregular points of Q with respect 
to dQ,i. 

As above, we deflne in f2 — Fj the minimal eikonal with one reflection in dQi such that 
Si^2{x,x) > Si^i{x,x). We construct an approximation 

2 / 

G{y, X, t) ~ Go{y, x,t) + J2Yl ^iAv^ t) 

k=l i=l 

where Gi^k{y,x,t) are ray solutions with eikonals Si^k{y,x) and Zi^k{y,x) chosen so as 
to minimize the boundary values of the sum at the boundary points of rays orthogonal 
to the boundary, as above. The trace of the double sum is calculated by introducing 
the change of variables 2rj fc(a;) = 5*4^^(0?, a;) and arclength Si{x) in dQi, as above. The 
Laplace evaluation of the integrals produces exponents that are the 2-periodic orbits in 

n. 

Eikonals with two or more reflections contribute exponents that are lengths of ex- 
tremal closed orbits with any number of reflections in the boundary, as in the case of 
simply connected domains. Thus the exponents (5„ in ()1.9|) consist of half the lengths of 
2-periodic orbits in Q and their multiples, and extremal lengths of closed periodic orbits 
with any number of reflections in the boundary and their multiples. 



5. Discussion 



First, we compare our result to that conjectured in 0. The leading exponent in the 
remainder of the short time expansion of P{t) in an asymptotic power series and that in 
the large s expansion of its Laplace transform are related by the well known formula ^Hj 



erfc 



2v^, 



where 



erfc (2;) 



for 2; ^ 1. 



Thus the exponential rate of blowup of the Laplace transform on a Stokes line is twice 
the square root of the exponential decay rate of the remainder term in the expansion 
of the trace. It follows from our result that the exponential blowup rate of the Laplace 
transform is twice the square root of the exponential rate of decay in the time domain. 

We illustrate our expansion for a disk, whose boundary has only one connected com- 
ponent and a single critical point. We consider points x = {xi,yi) and y = {x2,y2) inside 
a circle of radius R centered at the origin. The leading order eikonal is 



5*0 (y,x) = \x-y\ 
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When both x and y are on the x-axis, we have yi — y2 — ^ and Sq {y, x) — \xi — X2\- 
Denoting Xi = {xi, 0) and X2 = {x2, 0), we sec that the values of the eikonal on the x-axis 
arc So {xi, X2) — \xi — X2\. We assume that xi > 0. The boundary values of the eikonal 

are 

So {xi,X2) = R — xi at X2 = {R,0) 



anu 

So{xi,X2) = R + xi at CC2 = (— -R, 0). 
Thus the leading order ray approximation to Green's function G {y, x, t), 

Go{y,x,t)^2rt''^oiy,x)/4t^ 

misses the boundary conditions when x and y are on the x-axis, giving 

Go {x,, X2, t) = ^e-(^-^i)'/^* at X2 = {R, 0) 

and 

Go(a;i,a;2,t) = -^e-(^+^i)'/^* at 0^2 = (-i?, 0). (5.1) 

4:711 

The next eikonal, denoted {y, x), is given on the x-axis by 5*1 {xi, X2) = 2R—xi—X2, 
and its boundary values are 

Si{xi,X2) — R — xi at X2 — {R,0) 

Si {xi,X2) — SR — xi at X2 — {—R,0). 
approximation of Green's function G {y, x, t) , 

G {y, X, t) - Go {y, x, t) - Gi {y, x, t) , 

corresponding to the ray solutions Go (v, x, t) and 

Gi{x,y,t)^Zi{x,y,t)e-^i(y^^y'^t, 
will satisfy the boundary condition at {xi, R) if Zi [y, x, t) is chosen so that 

Zi{xi,X2,t) = ^ at CC2 = (-R,0). 

However, this approximation docs not satisfy the boundary condition at X2 — (— i?, 0). 
The error in the boundary values at X2 = {—R, 0) is 

Go {Xi, X2,t)-Gi {Xi, X2,t) = -Le-i^+-.r/^t_z^ ^^^^ ^^^^^ ^-4iR-.,f ^ ^ 

47rr 



and 

Thus the 
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and is of the same order of magnitude as that of the leading order approximation ()5.ip . 
To make up for the missed boundary condition the further approximation 



G (y, X, t) ~ Go (2/, X, t) - G, (y, x, t) - G^ {y, x, t) (5.2) 

can be used, with 

G,{y,x,t) = Z,{y,x,t)e-'"(y^'')/^^, 

where on the x-axis 

si (a?i, X2) = 2R + xi + X2 

and ^ 

Z2{xi,X2,t) = — at X2 = {-R,0). 

4:1X1 

This eikonal corresponds to rays with two reflections in the boundary. The approximation 
(15. 2|) decreases the error in the boundary condition at X2 = (— -R, 0) to 

-Z, {x,,X2,t)e-'^^-^^^'/'\ 
but misses the boundary condition at X2 = {R, 0) with error 

Goixi,X2,t)-Gi ixi,X2,t)-G2{xi,X2,t) = -Z2(a;i,X2,t)e-(^^+^'^)'/^* at X2 = (i?,0). 

This process gives successive approximations to Green's function with errors that decrease 
at transcendental rather than algebraic rates. 

The approximation to the trace produced by GQ{y,x,t) is the first algebraic term 
in the expansion p.7|) . The contributions of the terms —Gi {y,x,t) and —G2 {y,x,t) in 
the approximation ()5.2|1 of terms that are O (^y/ie~^^^^^ are identical, but with opposite 
signs and thus they cancel each other. The second term contributes a negative term 

that is O (v^e"^^"/*). The term O (v^e"^"/*) for small t corresponds to O (^-e~^^^^ 

for large positive s in the Laplace plane. The number 2R is the length of the periodic 
orbit of a billiard ball bouncing inside a circle with the center removed, that is, inside 
the domain — F, where the set of critical points F consists of the center. Similarly, 

the term O (^v^e~^^^/*j for small t corresponds to O ^-e~^^^^ for large positive s in the 

Laplace plane. The number 4R is the length of the minimal periodic orbit of a billiard ball 
bouncing inside a disk. We conclude that the conjecture of |H| should be supplemented 
with the orbit of length 2R. 

If Q is an annulus between two concentric circles, of radii a and b, respectively, {a > b), 
the two connected components of the boundary are the two circles and there are no critical 
points in the domain relative to either one of them. In this case 61 = {a — b). 

If Q is the ellipse 
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with a > b, the locus of critical points relative to the boundary is the segment 

a a 

on the X-axis. The segment F is the short diagonal of the evolute of the ellipse (the 
asteroid (ax)^^^ + {hy)^^^ = (a^ — b"^)"^^^). For the ellipse there are exponents in eq. p.9p 
which are 6i = 2b and its multiples and 62 = 2a and its multiples, as well as extremal 
periodic orbits with any number of reflections in the boundary. 

Finally, we observe that if the boundary is reflecting (i.e., a homogeneous Neumann 
boundary condition), the exponential decay rate of the transcendental terms in the expan- 
sion of the trace is the same as in the case of absorbing boundary (homogeneous Dirichlet 
boundary condition). In this case the second term in the expansion (|1.9p changes sign. 

Obviously, rays that are reflected from the boundary more than once also give rise to 
ray solutions. The number of ray solutions needed in the expansion ()3.7|) is determined 
by the required degree of asymptotic approximation of the boundary conditions. If only 
a finite sum of ray solutions satisfies the boundary conditions, the sum fl3.7|) is finite. 
Otherwise, additional ray solutions improve the degree of approximation of the boundary 
conditions, as described in the one- dimensional ray expansion in Section 2. 

Finally, the asymptotic convergence of the ray expansion follows from the maximum 
principle for the heat equation in a straightforward manner. 
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